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SUMMARY 


Bar specimens were cut from ingots of single crystal silicon, and acid- 
etched prior to testing. Artificial surface flaws were introduced in speci- 
mens by indentation with a Knoop hardness tester. The specimens were leaded 
in four-point bending to 95% of the nominal fracture stress, while keeping 
the surface area (Containing the flaw, wet with test liquids. Ho evidence of 
delayed fracture, and, therefore stress corrosion, of single crystal silicon 
was observed for liquid environments including water, acetone and aqueous 
solutions of NaCi, NH 4 OH, and HHO 3 , when tested with a flaw parallel to a 
(110) surface. The fracture toughness was calculated to be K^ c ■ 0.591 x 
10 6 N/m 3/2 . 
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DELAYED FRACTURE OF SILICON 


I. INTRODUCTION 

It is recognized that the fracture behavior of most Brittle Materials 
is markedly environment-sensitive; for example, delayed fracture (static 
fatigue) of a brittle element under imposed load (and residual stresses) 
may be caused by stress corrosion, or, by adsorption of an active species, 
and the resulting crack propagation. In the use of such elements in 
systems requiring significant structural performance, consideration of the 
effects of chemical environment on possible delayed fracture is of great 
importance. 

The purpose of this study was to investigate the effects of several en- 
vironments on delayed fracture of single crystal silicon. Silicon materials 
are receiving intensive study for use in photovoltaic energy conversion 
systems. The question of the sensitivity of silicon to stress corrosion 
with atmospheric water, organic liquids, acid and alkaline aqueous solu- 
tions, for instance, has not been previously- resolved, and may well be of 
considerable importance for the design of full-scale systems. In th„s in- 
vestigation bar specimens, cut from single crystal silicon, were to be sub- 
jected to static loading under several selected environments. 

II. PRIOR WORK 

Certain environment-sensitive fracture characteristics of ceramics glasses, 
and other brittle materials have been recognized for many years. It is 
well known that delayed fracture of silicate glasses and crystalline ceramics 
often is explained by stress corrosion reaction with atmospheric water since 
similar tests in vacuo or in completely dry air show an elimination or reduc- 
tion of delayed fracture. In studies of delayed fracture of an alumina 
ceramic, Chen and Knapp [l] found that the mean time to fracture, of stat- 
ically-loaded specimens, was Increased by about four orders of magnitude 
when the specimens were tested in dry argon rather than in moist air (50% 
relative humidity). An experimental procedure was adapted by Wiederhorn [2] 
to measure the propagation rate of cracks in glass and ceramics; he reported 
a thousand-fold increase in crack propagation velocity at a given stress in- 
tensity, when exposed to water instead of dry nitrogen. In studies of re- 
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lated properties, Westwood [3] has found a correlation between the hard- 
ness and zeta-potentlal of a non-metallic solid in a liquid environment, 
which indicates that the hardness is greatest when the zeta-potential is 
zero. Correspondingly, hardness is reported to be Influenced by the pH 
of an aqueous environment. Westwood also explained the variation of hard- 
ness of a soda-lime-silicate glass in aqueous solutions by the variation of 
induced surface potential. In addition, correlations are made [3] between 
the drilling rate and machining rates of non-metallic solids and the zeta- 
potential. The use of alcohols, in connection with the drilling or 
machining of ceramics and rocks , has been found advantageous and is ex- 
plained in terms of their effect on the zeta-potential. The effect of 
straigth chain organic alcohols and alkanes on crack propagation in various 
glasses also was studied by Freiman [4]. 

From the general aspects of environment-sensitive fracture, it can be 
assumed that surface-active environments will have little or no effect on 
the propagation of fast moving cracks in ceramic solids, but can be expected 
to influence greatly the propagation rate of slowly moving cracks. If the 
propagation rate is increased by an active environment, the effect usually 
is associated with (i) a stress corrosion process in which the stress en- 
hances the rate of corrosion at the crack tip relative to the sides, thus 
sharpening and deepening the crack, and/or with (ii) a lowering of the sur- 
face free energy by the adsorption of an active species, which reduces the 
energy contribution during fracture associated with new surface (Rebinder 
effect [5]). In the case of ceramics, the result of most prior studies 
have been explained best by assuming a stress corrosion process, but not 
all prior results can be understood in terms of this model. 

In the case of single crystal silicon, no prior investigations of delayed 
fracture are known to the authors. Jaccodine [8] has measured the fra.tured 
energy of silicon in the (111) planes, and has estimated values for the 
(100) and (110) planes. 

III. EXPERIMENTAL PROCEDURES 
Specimens 

The specimens used in this study were rectangular bars of single crystal 
silicon, approximately 2.54 cm x 0.30 cm x 0.20 cm in size. The specimens 
were cut, using a diamond saw (O.D. type), from ingots of single crystal 
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silicon, grown in-house at JPL by the Czochralski method, and mostly 
oriented as shown in Figure 1. A smaller number of specimens were tested 
with the (100) surface in tension. The surfaces of the specimens were treat- 
ed, prior to load-testing, by etching in an aqueous solution comprised of 
HN0 3 (15. 8N), HF (28. 8N), and HC 2 H 3 0 2 (17. 5N), in the proportions 1:2:3, 
for 1.0 minute at room temperature. In certain specimens, an artificial 
flaw was introduced in the mid-point of a tensile surface, following the 
method of Petrovic and Mendiratta [ 9 ], with a Knoop hardness-testing in- 
dentor. Figure 1 gives a schematic representation of the location and 
orientation of a crack produced by Knoop-indentation; when a lOOg load is 
used, the characteristic length of the major axis of the crack ranges from 
40-60yra. The indentation is always made in the surface to be under tensile 
stress, with the major axis of the crack normal to the direction of tensile 
stress. 

Load-testing 

(1) Short time loading 

Specimens both with and without indentiations were loaded in this part 
of the experimentation, using a four-point bending test fixture 
(Figure 2), with an Instron mechanical testing machine. The loading 
rate was 0.005'7min. Deflections under loading were recorded. 

(2) Delayed fracture 

Only indented specimens were used, for delayed fracture testing, with 
a modified four-point bending test fixture (Figure 3). A constant 
load, corresponding to 95% of the nominal fracture stress for short 
time loading, was applied for durations up to two weeks. During load- 
ing, the surface regions of the specimen containing the indentor-pro- 
duced crack was kept wet with a test liquid by means of a tissue wick, 
as shown in Figure 3. A timing clock, actuated by a contact micro- 
switch, monitored the time under load. 

(3) Incremental loading 

With certain specimens, a long period (more than one week) of static 
loading (at 95% of the nominal fracture stress) transpired without 
failure. To attempt to obtain failure data on incremental loading pro- 
cedure was initiated by a load producing 80% of the nominal fracture 

* Specimens were prepared and supplied by C.P. Chen, of the Jet 
Propulsion Laboratory. 
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Figure 2. The test fixture used for short 
time loadings. 




stress, followed by additional loads of 2.5% of the nominal fracture 
stress each 24 hours thereafter. This type of loading was used in an 
effort to cause delayed fracture by increasing the level of stress. 

Environmental conditions 

A a mentioned above, test liquids were applied to surface regions, con- 
taining artificial flaws, of specimens tested far delayed fracture, using 
wicks. Liquids used for this testing were: distilled water, aqueous so- 
lution of NaC&(10g/l), NH*0H(1.75M) , HNOj(2.76M), and acetone. All load- 
testing was conducted at room temperature. 


IV. EXPERIMENTAL RESULTS 

(1) Short time loading 

Fracture stress data, under short time loading, were obtained for 
specimens with artificial surface flaws (produced by Knoop inden- 
tation), and for specimens without artificial surface flaws. The 
results for ten specimens without indentations are gi,en in Table 
I and Figure 4, and indicate a mean modulus of rupture of 190.2 MN/m 2 
(27,5Sw psi), with a highest value of 281.9 MN/m 2 (40,890 psi), a 
lowest value of 133.8 MN/m 2 (19,400 psi), and a standard deviation of 
±40.1 MN/m 2 (5813 psi), or coefficient of variation of 21%. A typical 
fracture pattern for these specimens (without indentations) is shown 
in Figure 5. 

Fracture stress results for thirteen specimens with artificial surface 
flaws (Indented with a Knoop tester) are given in Table II and Figure 
6. Modulus of rupture values for indented specimens indicate a mean 
of 104.7 MN/m 2 (1^,190 psi), a high of 116.5 MN/m 2 (16,890 psi), a 
low of 89.4 MN/m 2 (12,970 psi) and a standard deviation of ±7.2 MN/m 2 
(1038 psi), or coefficient of variation of 6.8%. A typical fracture 
surface for indented specimens is shown in Figure 7. 

(2) Delayed fracture tests 

No significant evidence of delayed fracture was obtained for any of 
the specimens tested in this study. All specimens were indented before 
load-testing. Eleven specimens in contact with liquid water were 
loaded to 95% of the nominal (short time) fracture stress; although 
five of these specimen fractured immediately, upon application of the 
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TABLE I 


Fracture Stress of Silicon Specimens without 
Indentations under Short Time Loading 


Cumulative Fracture 


Rank 

Fracture 

Stress 

probability 

(i) 

(psi) 

(MN/m 2 ) 

(i/l+N) 

1 

19392 

133.70 

0.091 

2 

23787 

164.01 

0 182 

3 

24360 

167.96 

0.273 

4 

25166 

173.51 

0.364 

5 

25688 

177.11 

0.455 

6 

25945 

178.88 

0.545 

7 

2SJ32 

200.17 

0.636 

8 

30288 

208.83 

0.727 

9 

31265 

215.56 

0.818 

10 

40892 

281.94 

0.909 


190.18 


MN 

TD 7 


40.08 MN 




Mean 

Standard deviation 


(27,582 psi) 
( 5,813 psi) 





Figure 4. Fracture probability of specimens without Indentations 
under short time loading. 
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FIGURE 5. 


A typical fracture surface of a specimen, 
without an indentation, under short time 
loading. 



TABLE II 


Fracture Stress of Silicon Specimens with 
Indentations under Short Time Loading 


Cumulative Fracture 

Rank Fracture Stress probability 


(i) 

(psi) 


(MN/m 2 ) 

(i/l+N) 

1 

12969 


89.42 

0.071 

2 

13851 


95.50 

0.143 

3 

14751 


101.70 

0.214 

4 

14800 


102.04 

0.286 

5 

14834 


102.28 

0.357 

A 

14863 


102.48 

0.429 

7 

15290 


105.42 

0.5 

8 

15362 


105.92 

0.571 

9 

15668 


108.03 

0.643 

10 

15718 


108.37 

0.714 

11 

15980 


110.18 

0.786 

12 

16438 


113.34 

0.857 

13 

16894 


116.48 

0.929 


Mean = 

104.70 

(15,186 psi) 


Standard 

deviation = 

7.16 

^ ( 1,038 psi) 
m 2 
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Figure 6. Fracture probability of Indented specimens versus 
applied stress. 
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A typical fracture surface of a specimen, 
with an indentation, under a short time 
loading. 




load, the balance of the specimens continued to support the load 
after a period of two weeks. It was found that small vibrations, 
such as those produced by jarring the test stand, caused immediate 
fracture of some specimens. 

In addition, the several specimens tested each in contact with salt 
water, dilute solutions of NHi,OH and HNO3, and acetone, showed no 
indication of delayed fracture, and, therefore, no stress corrosion. 

(3) Incremental loading tests 

Results for the incremental loading of ten indented specimens in con- 
tact with water are given in Table III and Figure 6, and Indicate a 
mean modulus of rupture of 109.5 MN/m 2 (15,890 psi) , high and low 
values of 153.5 MN/m 2 (22,260 psi) and 85.8 MN/m 2 (12,450 psi) re- 
spectively, a standard deviation of ±17.9 MN/m 2 (2,600 psi), or co- 
efficient of variation of 16%. It may be noted that these results 
are similar to those obtained for short time loading of indented 
specimens (Table II). 

V. DISCUSSION OF THE RESULTS 

It was concluded, from these results, that no delayed fracture of single 
crystal silicon occurs, under an applied stress of 95% of the nominal 
fracture (short time) stress, for periods of two weeks in liquid environ- 
ments including water, acetone, and aqueous solutions of NaCJl, NH4OH and 
HNO3. Little, or no, stress corrosion appears to occur in such periods of 
study. It is interesting to compare the fracture probability curves 
(Figure 6) for short time loading, and for incremental loading (over ex- 
tended time periods) of specimens wetted with water. These curves are 
quite similar, and there is no evidence of a stress corrosion effect due 
to water. 

It is of additional interest to note certain fracture characterisitcs of 
single crystal silicon, especially in regards to the use of a Knoop in- 
dentor to introduce artificial cracks. The data for short time loading 
indicate, as may be expected, that the standard deviation of fracture 
stress is considerably reduced by introducing controlled artificial cracks 
with a Knoop indentor. The fracture surfaces of specimens show distinctive 
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TABLE III 


Frr.oture Stress of Silicon Specimens with Indendatlons 
under Incremental Loading with a Water Environment 


Rank 

Fracture 

Stress 

X of nominal 
fracture 

Cumulative Fracture 
probability 

(i) 

( 3i) 

(MN/m 2 ) 

% 

(i/l+N) 

1 

12445 

85.81 

82 

0.091 

2 

14248 

98.24 

94 

0.182 

3 

14792 

101.99 

97 

0.273 

4 

14964 

103.17 

99 

0.364 

5 

15247 

105.12 

100 

0.455 

6 

15309 

105.55 

100 

0.545 

7 

15571 

107.36 

103 

0.636 

8 

16780 

115.69 

110 

0.727 

9 

17251 

118.94 

114 

0.818 

10 

22259 

153.47 

147 

0.909 


Mean = 

MN 

109.33 ™ 
m 

(15,886 psi) 


Standard 

deviation = 

17.91 2? 
m 

( 2,597 psi) 





characteristics for specimens with, and without, indentations. A Knoop 
indentation produces a crack with a circular leading edge, as may be 
seen in Figure 8. 

The experimental data for short time load testing can be used to calcu- 
late the fracture toughness (K^) of single crystal silicon, using the 
following relation of Petrovic et al [10] : 

Kjc * o c MOra/Q)* 5 (D 

where a is the maximum outer-fiber tensile stress, M is a numerical 
c 

factor related to flaw and specimen geometry, a is the flaw depth, and 
Q is given by 

/ .w/2 

[sin 2 0 + cos 2 ©]* 5 d9 (2) 

o 

H 

where Q is the elliptic integral. For a typical indented specimen 
(specimen #14, Fig. 9) the values of 

Q = 1.6 

a * 22.1 x 10 ^m 

c ■ 21.0 x 10 ^m (major axis of the indentation) 
were obtained, leading to a calculation of the fracture toughness 

*10 * 0.591 x 10 6 N/m 3/2 

The Young's modulus of elasticity (E) for single crystal silicon was cal- 
culated, from the data for short time loadings, to be 

E * 8.96 x 10 4 MN/m 2 (13 x 10 6 psi) 

It is felt that the above value is lower than the true value, because the 
displacement under load was taken as the crosshead displacement of the 
testing machine. 

Examination of a number of fractured specimens confirmed that the plane of 
fracture was (111) in each case. Most specimens were loaded against the 
(110) surfaces, as shown in Figure 1QA, whereas certain specimens were 
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Figure 8. Schematic representation of a Knoop Indentation 
and the corresponding circular crack. 
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FIGURE 9. A semicircular crack produced by Rnoop 
indentation. 
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Figure 10. Schematic representation of the 

fracture planes (111) from loading 
on (A) (110), and (6) (100) surfaces. 



loaded against the (100) surfaces (Figure 10B); fracture followed (111) 
planes in both types of loading. 


VI. CONCLUSIONS 

It was concluded, from the results of this study, that 

(1) there is no evidence of delayed fracture (and stress corrosion) of 
single crystal silicon on liquid environments including water, ace- 
tone, and aqueous solutions of NaCfc, NHi»0H and HNO3 when tested with 
a flaw parallel to a (110) surface. 

(2) the fracture toughness of silicon, is calculated 

3/2 

to be K I( , - 0.591 x 10* N/m . 

(3) fracture occurs most readily along (111) planes, as shown by tests with 
specimens of different orientations. 

VII. RECOMMENDATIONS FOR FUTURE WORK 

(1) Although no evidence of delayed fracture, and, therefore stress corro- 
sion, was observed for the environmental influences of this study, it 
seems likely that corrosive liquids for silicon remain to be identi- 
fied. Additional liquids, such as other organic liquids, acids and 
alkalis should be evaluated. 

(2) The fracture of silicon under several stress states, including com- 
pression, and combined stress, should be studied. Possible types of 
loading may include that producing twisting, and one developing 
equlbiaxial tension. 

(3) The influence of flaw orientation on the fracture of silicon is a 
study topic of importance. 

(A) The fracture of silicon over a range of temperatures should be investi- 
gated. 
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